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1 Apologia 

In these informal Lecture Notes on formalism, I gather together a few elements of covariant linear 
perturbation theory and its applications to inflation, dark energy or matter and the cosmic mi- 
crowave background (CMB). I make no attempt to cite properly the original source material. My 
usual, lame defense is "it's all in Bardeen's paper PP" - go read it and ignore these lecture notes! 

Seriously though, in addition to PP, I have heavily relied on a few sources in compiling the 
various sections and refer the reader to references therein. The covariant formalism presented in ^ 
and gauge transformation material in ^ draws from [2] distilled in Applications to inflation in 
^ draw from |lj , to dark energy and matter in ^ from , and to the CMB in ^ from fS] . For the 
last topic, despite great recent progress on the phenomenology, I have limited myself to the formal 
aspects that relate to covariant perturbation theory. My defense here is "everything that doesn't 
go out of date as soon as it's written is in Peebles & Yu [7j" - go read it and stop looking for CMB 
reviews! My other defense is Matias Zaldarriaga was supposed to cover that in this School so blame 
him! As for the equation density and opaqueness of these Lecture Notes, I have no excuse - take 
it as an homage to the Dick Bond lectures of my own student days - at least there are no figures 
with 100 curves. As Martin White would say, you are physicists; you don't need figures! 

2 Formalism 

2.1 Covariant Approach 

Perturbation theory proceeds by linearizing the Einstein equations 

G^, = SttGT^,, (1) 

around a background metric. Here Gf^i, is the Einstein tensor and T^u is the stress energy tensor. 
The Bianchi identity V ^G^'^ = guarantees the covariant conservation of the total stress-energy 

V^T'^^ = 0. (2) 

The conservation equation is redundant with the Einstein equation but is a particularly useful rep- 
resentation of the equations when there are multiple components to the matter that are separately 
conserved. 

Let us begin by distinguishing between covariant equations and invariant variables: 

Covariant = equations takes same form in all coordinate systems 
Invariant = variable takes the same value in all coordinate systems 
In an evolving universe, the meaning of a density perturbation necessarily requires the specification 
of the time slicing in relation to the background and hence there is no such thing as a gauge 
or coordinate invariant density perturbation. Likewise for elements of the metric. The general 
covariance of the Einstein and conservation equations, on the other hand, can be preserved. With 
covariant equations one can, after the fact of their derivation, choose the gauge that best suits a 
given physical problem, e.g. the evolution of inflationary, dark energy, or CMB fluctuations. 

2.2 Metric Representation 

Let us take the background metric to be the general homogeneous and isotropic, or Friedmann- 
Robertson- Walker form. Here the degrees of freedom are the (comoving) spatial curvature K and 
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an overall scale factor for the expansion a such that the line element 

ds^ = g^udx^dx"" = a^{-dif + ^ijdx'dx^) . (3) 

The time variable is 77 the conformal time and we normalize the scale factor to a = 1 today, 
r]{a = 1) = r/o- The three metric ^ij can be represented in spherical coordinates as 

-fijdx'dx^ = dD^ + D\dQ , (4) 

where D is the comoving distance and Da = K~^/'^ sin(iC^/^D) is the angular diameter distance. 
A general perturbation to the FRW metric may be represented as 

g'' = -a-\l-2A), 

g^^ = a-2(y^- - 2Hlj'^ - 2Hj^) , (5) 

yielding (1) A = a scalar potential; (3) a vector shift; (1) Hi a scalar perturbation to the spatial 
curvature; (5) a trace-free distortion to spatial metric, for a total of (10) degrees of freedom. 
This gives a complete representation of the symmetric 4x4 metric tensor. 



2.3 Matter Representation 

Likewise expand the matter stress energy tensor around a homogeneous density p and pressure p 
by introducing 10 degrees of freedom 

T% = -p-5p, 

A = {p + p){v^-Bi), 

V = -{p + py, 

T) = {p + 6p)6'^+pll'^, (6) 

yielding (1) 5p a scalar density perturbation; (3) Vi a vector velocity, (1) Sp a scalar pressure 

perturbation; (5) Hij a tensor anisotropic stress perturbation. So far the treatment of the matter 
and metric is fully general and applies to any type of matter or coordinate choice including non- 
linearities in the matter, e.g. cosmological defects. 



2.4 Closure 

The Einstein equations do not specify a closed system of equations and require the addition of 
supplemental conditions provided by the microphysics to fix the dynamical degrees of freedom. In 
the simplest cases this closure is established through equations of state for the background and 
perturbations. 

In counting the dynamical degrees of freedom, one must recall that the conservation equations 
are redundant with the Bianchi identities and that 4 degrees of freedom are absorbed by the choice 
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of coordinate system 

20 Variables (10 metric; 10 matter) 
— 10 Einstein equations 
—4 Conservation equations 
+4 Bianchi identities 

Gauge (coordinate choice 1 time, 3 space) 



-4 



6 Degrees of freedom 

Without loss of generality these six dynamical degrees of freedom can be taken to be the 6 compo- 
nents of the matter stress tensor. 

Since the background is isotropic, the unperturbed matter is described entirely by the pressure 
p{a) or equivalently 'w{a) = p{a)/p{a), the equation of state parameter. For the perturbations 
specification of the relationship between Sp, H and the other perturbations, e.g. Sp and v suffices. 

2.5 Friedmann Equations 

The unperturbed Einstein equation yields the ordinary Friedmann equation which relates the ex- 
pansion rate to the energy density 

where overdots represent conformal time derivatives, and the acceleration Friedmann equation 
which relates the change in the expansion rate to the equation of state 

d /d\ AttG 2, 



,„\a,' + (8) 

SO that w = p/ p < —1/3 implies acceleration of the expansion. 
The conservation law VT^j^ = implies 

- = -3(1 + «;)-, (9) 
p a 

which by virtue of the Bianchi identity can be derived from the Friedmann equation. 
The counting exercise for the background becomes 

20 Variables (10 metric; 10 matter) 

—17 Homogeneity and Isotropy 

—2 Einstein equations 

—1 Conservation equations 

-1-1 Bianchi identities 



1 Degree of freedom 

Without loss of generality this degree of freedom can be chosen to be the equation of state w{a) = 
p{a)/p{a). 
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2.6 Linecirization and Eigenmodes 

For the inhomogenous universe, the Einstein tensor Gfj,^, is in general constructed out of a nonhnear 
combination of metric fluctuations. Metric fluctuations typically remain small even in the presence 
of larger matter fluctuations. Hence we linearize the left hand side of the Einstein equations to 
obtain a set of partial differential equations that are linear in the variables. 

These equations may then be decoupled into a set of ordinary differential equations by employing 
normal modes under translation and rotation. The scalar, vector and tensor eigenmodes of the 
Laplacian operator form a complete set 

V2Q(0) = -fc2Q(0) S, 

v2Q,(±^) = -eql^'^ V, (10) 

V^Qjf) = -PQlf T. 
In a spatially flat {K = 0) universe, the eigenmodes are essentially plane waves 

= exp(ik-x), 
Qi^'^^ = ^^2)i exp(ik • x) , 

/ 3 

Q\f^^ = -y-(ei±ie2)i(ei±ie2)jexp(ik-x), (11) 

where ei , 62 are unit vectors spanning the plane transverse to k. 

Vector modes represent divergence- free vectors (vorticity); tensor modes represent transverse 
traceless tensors (gravitational waves) 

V'q[^^^=0, V'Qlf)=0, 7*^g(f)=0. (12) 

Curl free vectors and the longitudinal components of tensors are represented with covariant deriva- 
tives of the scalar and vector modes 

Qr = -^[V.Qf^+V,QS^^)], (13) 

For the fcth eigenmode, the scalar components become 

^(x) = A{k)Q^^), Hl{^) = HL{k)Q(^\ 

the vectors components become 

1 1 
B,(x) = ^S(-)(A;)Qf), Vi{yL) = J] i;^"*) (fc) ^ , (15) 

m=— 1 m=— 1 

and the tensors components become 

Ht^M)= E , n,,(x)= nM(fc)Q;7). (16) 

m=— 2 m=— 2 

An arbitrary set of spatial perturbations can be formed through a superposition of the eigenmodes 
given their completeness. 



(14) 
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2.7 Coveiriant Scalar Equations 

The Einstein equations for the scalar modes (suppressing superscripts) become 



6 a K 



fc2 



)] 



47rGa2 
-87rGa^pn , 



5p + ^-{p + p)- 



B 



B 



a 
2- 
a 



2 - 



a d 
a dr) 



3 





d a 


A- 


dr) a 



4TrGa^{6p+ -5p) 



The conservation equations become the continuity and Navier-Stokes equations 



d d 
dr] a 



d , d 
+ 4_ 

dr] a 



bp + 2)— bp 
a 

B 



ip + p)- 



k 



~{p + p){kv + 3HL), 
bp-^il-3^)pU + {p + p)A. 



(17) 

(18) 
(19) 



In the absence of anisotropic stress, the Navier-Stokes equation is known as the Euler equation. 

These equations are not independent due to the Bianchi identity. For numerical solutions it 
suffices to retain 2 Einstein equations and the 2 conservation laws. The counting exercise becomes 

8 Variables (4 metric; 4 matter) 
—4 Einstein equations 
—2 Conservation equations 
+2 Bianchi identities 

—2 Gauge (coordinate choice 1 time, 1 space) 
2 Degrees of freedom 

Without loss of generality we may choose the dynamical components to be those of the stress 
tensor bp, U. Microphysics defines their relationship to the density and velocity perturbations. In 
fluid dynamics this is the familiar need of a prescription for viscosity and heat conduction (entropy 
generation) to close the system of equations. 



2.8 Coveiriant Vector Equations 

The vector Einstein equations become 



# + 2- 

dr] a 

and the conservation equations become 



(fc5(±l)-ij(,±^)) 



167rGa'^{p + p)- 



/±i) _ ^(±1) 



k 



-87rGa2pn(±^ , 



d a 
— +4- 
dr] a 



iP + P)- 



/±i) _ ^(±1) 
k 



= 4(i-4)^n(-). 



(20) 



(21) 
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Since gravity provides no source to vorticity, any initial vector perturbation will simply decay unless 
vector anisotropic stresses are continuously generated in the matter. In the absence of nonlinearities 
in the matter, e.g. from cosmological defects, vector perturbations can generally be ignored. 
The counting exercise for vectors becomes 

8 Variables (4 metric; 4 matter) 

—4 Einstein equations 

—2 Conservation equations 

+2 Bianchi identities 

—2 Gauge (coordinate choice 2 space)) 



2 Degrees of freedom 

Without loss of generality, we can choose these to be the vector components of the stress tensor 

n(±i). 

2.9 CovEiriant Tensor Equations 

The Einstein equation for the tensor modes is 



A 2+2-:^ + (/c' + 2i^) 



4^') = 87rGa2pn(±2) , (22) 



and the counting exercise becomes 



4 Variables (2 metric; 2 matter) 

—2 Einstein equations 

—0 Conservation equations 

+0 Bianchi identities 

—0 Gauge (coordinate choice 1 time, 1 space) 



2 Degrees of freedom 

without loss of generality we can choose these to be represented by the tensor components of the 

stress tensor 

In the absence of anisotropic stresses and spatial curvature, the tensor equation becomes a 
simple source-free gravitational wave propagation equation 

4^^)+2^4±^)+/c24±2)=0, (23) 

which has solutions 

H^^^\kr]) = CiHiikf]) + C2H2{kri) , (24) 
Hi{x) oc x"™im(x) , 
H2{x) oc x~'^nm{x) , 
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where m = (1 — 3w)/{l + Sw). If w > —1/3 then the gravitational wave amphtude is constant 
above horizon x <C 1 and then oscillates and damps. If w < —1/3 then gravity wave oscillates 
and freezes into some value. It will be useful to recall these solutions when considering the Klein 
Gordon equation for scalar field fluctuations during inflation and dark energy domination. 

2.10 Multicomponent universe 

With multiple matter components, the Einstein equations of course remain valid but with the 
associations 

J 

p^M=^pX™^ (25) 

J 

where J indexes the different components, e.g. photons, baryons, neutrinos, dark matter, dark 
energy, inflaton, cosmological defects. The conservation equations remain valid for each non- 
interacting subsystem. Interactions can be represented as non-ideality in the specification of the 
stress degrees of freedom (e.g. viscous and entropic terms) in each subsystem or by explicitly 
writing down energy and momentum exchange terms in separate conservation equations (see e.g. 



3 Gauge 
3.1 Semantics 

The covariant equations of the last section hold true under any coordinate choice for the relation- 
ship between the unperturbed FRW background and the perturbations. Choice of a particular 
relationship, called a gauge choice can help simplify the equations for the physical conditions at 
hand. The price to pay is that the perturbation variables are quantities that take on the meaning 
of say a density perturbation or a spatial curvature perturbation only on a specific gauge - but 
that is in any case unavoidable. 

Since the preferred gauge for simplifying the physics can change as the universe evolves, it 
is often useful to access the perturbation variables of one gauge from those of another gauge. 
This operation proceeds by writing down a covariant form for gauge specific variables through 
the properties of gauge transformation. It is the analogue of deriving covariant equations for the 
dynamics but for the perturbation variables themselves. Since these gauge specific variables are 
now thought of in a coordinate independent way, this procedure is called in the literature endowing 
a variable with a "gauge invariant" meaning. Note neither the numerical values nor the physical 
interpretation of the variables has changed. "Gauge invariance" here is an operational distinction 
that indicates a freedom to calculate gauge specific quantities from an arbitrary gauge. Perturbation 
variables still only take on their given meaning in the given gauge. Under this definition of gauge 
invariance, the perturbation variables of any fully specified gauge is gauge invariant. We will 
hereafter avoid using this terminology. 
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To reduce the opportunity for confusion, we will name the metric perturbations in the various 
common gauges separately. We will however rely on context to distinguish between matter variables 
on the various gauges. 



3.2 Gauge Transformation 

In an evolving inhomogcneous universe, metric and matter fluctuations take on different values in 
different coordinate systems. Consider a general coordinate transformation 

ri = fi + T, (26) 

x' = x' + L\ 

Under this general transformation, the metric and stress energy tensors transform as tensors. The 
elements of these tensors are the perturbation variables and their numerical values change with the 
transformation. 

The coordinate choice represented by (T, U) can be decomposed into scalar and vector modes. 
With L^^^ = L, the scalar mode variables transform as 

A = A-f--T, 
a 

B = B + L + kT, 

Hl = Hl — -L — -T , 
6 a 

HT = HT + kL, (27) 



for the metric and 



Spj = Spj - PjT, 

Spj = dp J - pjT, 
vj = vj + L (28) 



for the matter. 

For the vector mode variables 

^(±1) ^^(±i)+^(±i)^ 

4±i)=4±i)+L(±i). (29) 

The tensor mode variables are invariant under the gauge transformation as are the components of 
the anisotropic stress tensor. 

A gauge is fully specified if there is an explicit prescription for (T, L) for scalars and {L^^^^^ ) for 
the vectors to get to the desired frame. Gauges are typically defined by conditions on the metric 
or matter fluctuations. We now consider several common scalar gauge choices and their uses. 
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3.3 Newtonian Gauge 

The Newtonian or longitudinal gauge is defined by diagonal metric fluctuations 

B = Ht = 0, 

= A (Newtonian potential) , 
$ = Hl (Newtonian curvature) , 



(30) 



where ^ plays the role of the gravitational potential in the Newtonian approximation and ^ is 
the Newtonian spatial curvature. This condition completely fixes the gauge by giving explicit 
expressions for the gauge transformation 



L 
T 



B I d ~ 



(31) 



The Einstein equations become 



/ \V 



(32) 



Note that for scales inside the Hubble length k{a/a)^^ » 1, the first equation becomes the Poisson 
equation for the perturbations and if the anisotropic stress of the matter is much less than the 
energy density perturbation (as in the case of non-relativistic matter), ^ k, — $. 
The conservation laws for the Jth non-interacting subsystem becomes 



d d 

drj a 

d , d 
+ 4_ 

dr] a 



3-Spj 
a 



(PJ+PJ) 



vj_ 
k 



-{pj+Pj){kvj + 3^), 
2 K 

^pj - 3(1 - 3p)wnj + {pj+pj)-^ . 



(33) 
(34) 



Aside from the ^ term, these are the usual relativistic fluid equations. Since ^ represents a 
perturbation to the scale factor, $ represents the perturbation to the redshifting of the density in 
an expanding universe [see Eqn. Q]. 

The Newtonian gauge is useful since it most closely corresponds to Newtonian gravity. A 
drawback is that as k{d/a)^^ — > there are relativistic corrections to the Poisson equation which 
make a straightforward implementation of the equations numerically unstable. Likewise in this 
limit, the metric effects on density though <I> muddle the interpretation of the conservation law. 

3.4 Comoving Gauge 

The comoving gauge is defined by 

B = v iT^ = 0), 
Ht = 0, 
C = A, 

C = Hi (comoving curvature) , (35) 
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which completely fixes the gauge through 



The Einstein equations become 



T={v- B)/k , 
L = -Hr/k . 



(36) 



+ 2-v + k{C + = -STrGa^pn , 



and the conservation laws become 



df] a 



Spj + 3-Spj = -{pj+ pj){kvj + 3C) 
a 



d d 
+ 4_ 

dri a 



(PJ+PJ) 



k 



Spj - ^(1 - 3^)pjnj + (pj+pj)^. 



(37) 

(38) 
(39) 



In particular the Navier-Stokes equation for the total matter becomes an algebraic relation between 
total stress fluctuations and the potential 



(p + p)C = -^P+o 1 



3K\ 
TJ 



pU 



(40) 



so that these equations are a complete set. 

Eliminating ^ allows us to write a conservation law for the comoving curvature 



C + Kv/k = - 



a \_ p + p 3 \ k"^ 



p 



-n 
p+p . 



(41) 



On scales well below the curvature scale, this equation states that the comoving curvature only 
changes in response to stress gradients which move matter around. This statement corresponds to 
the non-rclativistic intuition that causality should prohibit evolution above the horizon scale. This 
conservation law is the fundamental virtue of the comoving gauge. An auxiliary consideration is 
that the comoving gauge variables are numerically stable. 

Prom an alternate gauge choice, one can construct the Newtonian variables via the gauge 
transformation relation 



6 a k 



B 



(42) 



For example, since in the Newtonian gauge Ht and v take on the same values 

a k 

:[*-*/(«/«)]■ 



Sl + wl + K{d/a)-^' 



(43) 



If the curvature is negligible, stress perturbations are negligible, and the equation of state w is 
constant 



3 + 3w 
5 + 3w 



C: 



(44) 
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so that <I> tends to be of the same order as but changes when the equation of state changes. 
The Newtonian curvature can also be obtained through the comoving density perturbations via 



(45) 



comoving 

which has the added benefit of taking the form of a simple non-relativistic Poisson equation. This 
Poisson equation allows us to rewrite the conservation equation for the comoving curvature in the 
case of negligible background curvature and anisotropic stress as 

= k- (a/a) — 

5p 



din a 



comoving 



5 + 3u; ^ ' ' Sp 



comovmg 



(46) 



For adiabatic stresses where Sp = c^Sp, the change in the comoving curvature is negligible for 
Cak{a/a)''^ ~ Cakr] <C 1, i.e. on scales below the sound horizon. 

3.5 Synchronous Gauge 

The synchronous gauge confines the metric perturbations to the spatial degrees of freedom 

A = B = Q, 



a ^ I dr]aA + cia ^ 



J dri{B 



+ kT) + C2 , 



(47) 



The metric conditions do not fully specify the gauge and need to be supplemented by a definition 
of (ci,C2). Usually one defines Ci through the condition that the dark matter has zero velocity 
in the initial conditions and C2 through the setting of the initial curvature perturbation. With a 
completely specified gauge condition, the synchronous gauge is as valid as any other gauge. The 
variables rjr and hz form a stable system for numerical solutions and hence the S3mchronous gauge 
has been extensively used in numerical solutions. 
The Einstein equations give 

K 

VT - 



■• a ■ 

hi + -hi 
a 



-87rGa2(5p + 3(5p) , 



while the conservation equations give 



d a 
dr] a 
' d_ 
drj 



Spj + = -{PJ +Pj){kvj + ^/il) , 



4^ 



2 K 
Spj ~ 3(1 ~ ^]^)PJ^J- 



(48) 

(49) 
(50) 



Note that the lack of a potential in synchronous gauge implies that there are no gravitational forces 
in the Navier-Stokes equation. Hence for stress free matter like cold dark matter, zero velocity 
initially implies zero velocity always. 
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3.6 Spatially Flat Gauge: 

Conversely the spatially flat gauge eliminates spatial metric perturbations 

Hl = Ht = 0, 
ap = A , 

Pf = B, 



T 
L 



, -1 

a 
a 

Hr/k, 



Hl + -Ht 



The Einstein equations give 



a 



/3f + 2-[3f + kap = -8TrGa^pU/k , 
a 

-ap — rPF = 47rGa [p + p) — ; 

a k k 



and the conservation equations give 



d a 
ar] a 



d , d 
+ 4_ 

dr] a 



(PJ + PJ 



5pj + 3-5pj 
a 

vj- f^F 



(51) 



k 



-{pj +Pj)kvj , 
2 K 

^PJ ~ ~ ^Ja^PJ^J + {PJ +Pj)<^F ■ 



(52) 

(53) 
(54) 



The spatially flat gauge is useful in that it is the complement of the comoving gauge. In particular 
the comoving curvature is constructed from Eqn. (|42() as 

dv — 13f 



c 



(55) 



a k 

This gauge is most often used in inflationary calculations where v — Pf is closely related to pertur 
bations in the inflaton field. 



3.7 Uniform Density Gauge: 

Finally, one can eliminate the density perturbation 6p = with the choice 

Ht = 0, 
Cs 
Bs 
As 

T 



The Einstein equations become 



L 



ik' - 3K)[Cs + 

a k 

-^5 - Cc -r^& 

a k 



B 
A 
6p 
P 

-Hx/k 



v-Bs 



UnGa^-ip + p) , 
a k 

Bs 



A-KGa^ip + p)- 



(56) 



(57) 
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from which As may be ehminated in favor of Q and Bg. The conservation equations become 

" d 



a 

ar] a 



d , d 
+ 4_ 

dr] a 



{PJ + PJ 



6pj + 3-5pj 
a 

vj - Bs 



-{pj +Pj){kvj + 3(5) , 
Spj - ^(1 - 3^)pjnj + ipj+pj)As . 



(58) 



k 3 

Notice that the continuity equation for the net density perturbation becomes a conservation law 

d 6p 



Cs 



a p + p 



—kv . 
3 



(59) 



Furthermore since 5p = 0, 5p is the non-adiabatic stress [see Eqn. (|46|) ] and the conservation law 
resembles that of the comoving curvature [see Eqn. (|41|) ]. More specifically, the two curvatures are 
related by 

1 5p 



C<5 = C + T7: 



comovmg 



(60) 



3(P + P) 

By the same argument as that following Eqn. H46|) . these two curvatures coincide outside the horizon 
if the stresses are adiabatic. Hence they are often used interchangeably in the literature. It bears 
an even simpler relationship to density fluctuations in the spatially flat gauge 



1 5p 



3 (p + p) 



fiat 



(61) 



For a single particle species 5p/{p + p) = 5n/n^ the number density fluctuation. These simple 
relationships make Q and related variables useful for the consideration of relative number density 
or isocurvature perturbations in a multicomponent system. 



4 Inflationary Perturbations 

4.1 Horizon, Flatness, Relics Redux 

Inflation was originally proposed to solve the horizon, flatness and relic problem: that the cosmic 
microwave background (CMB) temperature is isotropic across scales larger than Hubble length 
at recombination; that the spatial curvature scale is at least comparable to the current Hubble 
length; that the energy density is not dominated by defect relics from phase transitions in the early 
universe. 

Measurements of fluctuations in the CMB imply a stronger version of the problems. They 
imply spatial curvature fluctuations above the Hubble length and rule out cosmological defects as 
the primary source of structure in the universe, i.e. they are absent as a signiflcant component of 
the perturbations as well as the background. The conservation of the comoving curvature above 
the comoving Hubble length (d/a)^^ <^ 1 leaves a paradox as to their origin. As long as the Hubble 
length monotonically increases, there will always be a time r] ~ A;"^ before which C, cannot have 
been generated dynamically. For sufficiently small A;, this time exceeds the recombination epoch at 
a ~ 10~^ where the CMB fluctuations are formed. 

The comoving Hubble length decreases if p scales more slowly than a^^ 01 w < —1/3, i.e. if 
the expansion accelerates [see Eqn. (|7l8p ]. Vacuum energy or a cosmological constant can provide 
for acceleration but not for re-entry into the matter-radiation dominated expansion. On the other 
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hand, a scalar field has both a potential and kinetic energy and provides a mechanism by which 
acceleration can end. A major success of the inflationary model is that the same mechanism that 
solves the horizon, flatness and relic problems by making the observed universe today much smaller 
than the Hubble length at the beginning of inflation predicts super-Hubble curvature perturbations 
from quantum fluctuations in the scalar field. The evolution of scalar field fluctuations can be 
usefully described in the gauge covariant framework. 

4.2 Scalar Fields 

The covariant perturbation theory described in the previous sections also applies to the inflationary 
universe. The stress-energy tensor of a scalar field rolling in the potential is 

n = V'^<^V,¥'-^(V>Va(^ + 2F)^'',. (62) 

For the backround (p = and 

P^ = \a-m + V, V4>=\a~m-V. (63) 

If the scalar field is kinetic energy dominated tf^ = p,f,/p,), 1, whereas if it is potential energy 
dominated =P4>/ P4> ~1- 

The equation of motion for the scalar field is simply the energy conservation equation 

P0 = -3(p</, +p</,)^ , (64) 



re-written in terms of and V 
Likewise for the perturbations = 0o + 



</)o + 2-(^o + aV = 0. (65) 
a 



--a-\^^^i-4HiA) + V'4>i, 
--a-\Mi-4^lA)-V'4>^, 

P(l)7f<j> = , (66) 



the continuity equation implies 



6i = -2-<j)i - (k^ + aV")<^i + iA- 3Hl - kBUo - 2Ao?V' , (67) 
a 



while the Euler equation expresses an identity. 



4.3 Gauge Choice 

We are interested in inflation as a means of generating comoving curvature perturbations and so 
naively it would seem that inflationary perturbations should be calculated in the comoving gauge. 
However for a scalar field dominated universe, the comoving gauge sets B = v = and hence 
01 = 0. In the comoving gauge, the scalar field carries no perturbations by definition! This fact 
will be useful for treating scalar fields as a dark energy candidate in the next section: in the absence 
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of scalar field fluctuations, energy density and pressure perturbations come purely from the kinetic 
terms so that Spif, = 5p(j, yielding stable perturbations within the Hubble length. 

The gauge covariant approach allows us to compute the comoving curvature from variables in 
another gauge. Since a scalar field transforms as a scalar field 



'>! = <t>i- (t>oT 



(68) 



the comoving gauge is obtained from an arbitrary gauge by the time slicing change T = ^i/^q. 
The comoving curvature becomes 



6 a 



= Hl + 



Hx a (pi 



(69) 



and hence the simplest gauge from which to calculate the comoving curvature is one in which 
ifi = Ht = 0, i.e. the spatially flat gauge. In this case 



(70) 



and so a calculation of 4>i trivially gives (. Notice that the proportionality constant involves ^q. 
The slower the background field is rolling, the larger the curvature fluctuation implied by a given 
field fluctuation. The reason is that the time surfaces must be warped by a correspondingly larger 
amount to compensate cpi. We will hereafter drop the tildes and assume that the scalar field 
fluctuation applies to the spatially flat gauge. 

4.4 Perturbation Evolution 

In the spatially fiat gauge the scalar field equation of motion can be written in a surprisingly 
compact form. Beginning with the spatially fiat gauge equation 



.a ■ 



(f)^ = -2-<pi - (k^ + a^V'Ui + (ap - kppUo - lapa^V 
a 

the metric terms may be eliminated through Einstein equations 

-1 



ap = AttGo^ ( - ) (/90 + Pcj,){v^ - PF)/k 



and > \K\ 



A;/?F = 47rGa2 ( - 
= 47rG 



.0, 



a 

^) {<j>o4>i + aV>i) - {AttG^o?Mi + 30001 



(71) 
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so that ap, ap — kfSp oc with proportionaUty that depends only on the background evolution, 
i.e. the Einstein and scalar field equations reduce to a single second order differential equation with 
the form of an expansion damped oscillator 

<Ai + 2-,/.i + [k" + . (72) 

a 

The expression for f[r}) can be given explicitly in terms of the parameters 

e^-{l + w^)= (73) 
2 1 + ^(pf^/a^V 

which represents the deviation from a de Sitter expansion and 

■ \ -1 

a) P") 

which represents the deviation from the overdamped limit of (P(f)Q/dt^ = 0, where dt = adrj. When 
small, these quantities are known as the slow-roll parameters. The Friedmann equations become 



2 

.2,-1 



d a\ a 



(1 - 6) , (75) 



drj \a J \a 
and the background field equation becomes 

,/.o-(3 + <5) = -aV, (76) 
a 

Together they yield the equation of motion for e 

e = 2e(J + e)-. (77) 
a 

It is straightforward now to show that defining u = a(j) the perturbation equation becomes 

u+[k^ + gi7])]u = 0, (78) 

where 

9{^) ^ f{ri) + e - 2 = - f - y [2 + 35 + 2e + (5 + 6)(5 + 26)] - -5 

\a a 



(79) 



and ^ 

^ = a 4>o- (80) 

For any background field, Eqn. (|78jl can be solved to yield the evolution of the scalar field fiuctua- 
tion. 
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4.5 Slow Roll Limit 

In the slow roll limit e, 5 <C 1 and the calculation simplifies dramatically. The slow roll parameters 
are usually written in terms of the derivatives of the potential 

m/a^V 1 (V 



(81) 



1 + Ul/a?V 167rG V V 



. 1 V" 



(^0 Va / SvrG V 

The slow roll limit requires a very flat potential. 

In the slow roll limit, the perturbation equation is simply 

2 



n + [r - 2 { - ) ]u = . (82) 



With the conformal time measured from the end of inflation 

fi = r]- ??end , 



where aH = a/a, the equation becomes 



il+lk^ - = . (84) 



This equation has the exact solution 



u = A{k + -)e-'^^ , (85) 

where ^ is a constant. For \kfj\ 3> 1 (early times, inside Hubble length) the field behaves as free 
oscillator 

lim u = Ake-'^^. (86) 

Ifcj)!— >oo 

For \kf}\ <C 1 (late times, 3> Hubble length), the fiuctuation freezes in 

lim u = —A = iHaA , 
^1 = iHA , 

'a\ ^_ 

00 

which in the slow-roll approximation can be simplified through 

1 SirGa^V 3 AttG 4tt 



C = -iHA 



(87) 



yielding the comoving curvature power spectrum 

,,^fc3|^|2 2fc3 

^ 27r2 TT empj 

All that remains is to set the normalization constant A through quantum fiuctuations of the free 
oscillator. 
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4.6 Quantum Fluctuations 

Inside the Hubble length, the classical equation of motion for u is the simple harmonic oscillator 
equation 

u + k'^u = 0, (89) 

which can be quantized as 

u = u{k,ri)a + u*{k,ri)a^ , (90) 
and normalized to zero point fluctuations in the Minkowski vacuum [u,du/dr]] = i, 



Thus A = (2A;^)^/^ and curvature power spectrum becomes 

1 



Al^-^. (92) 

The curvature power spectrum is scale invariant to the extent that H is constant during inflation. 
Evolution in H produces a tilt in the spectrum 

dlnA^ _ 
din A; ^ 

dlnA; dlnA;' ^ ^ 

evaluated at Hubble crossing when the fluctuation freezes 

dlnH. kdH. dfj 



(94) 

= ^(--Hh)\l^=-e, (95) 



d In A; ' -^=^=1 H dfj ' -''^=^ dk ' -^^=i 



where aH = —1/fj = k. Finally with 



dine dine ^.^ ^.^ , 

—— = --— = -2{6 + e)-v = 2{S + e), 96 
dlnA; dln77 a 

the tilt becomes 

ns = 1 - 4e - 2(5 (97) 

in the slow-roll approximation. 
4.7 Gravitational Waves 

Any nearly masslcss degree of freedom will acquire quantum fluctuations during inflation. The 
inflaton is only special in that it carries the energy density of the universe. Other degrees of 
freedom result in isocurvature perturbations. In particular consider the gravitational wave degrees 
of freedom. Their classical equation of motion resembles the scalar field equation 

4^^)+2^4±^)+fc24±^)=0, (98) 
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and hence acquires fluctuations in same manner as (f). Setting the normahzation with 

the gravitational wave power spectrum in each component (polarization components Hj, = (/i+± 
ihx)/VQ) is 

The gravitational wave fluctuations are scale invariant with power oc i?^ oc V oc Ef where Ei 
is the energy scale of inflation. A detection of gravitational waves from inflation would yield a 
measurement of the energy scale of inflation. 

Finally the tensor tilt 

dlnA?. dlnH „ 

which yields a consistency relation between tensor-scalar ratio and tensor tilt 

5 Dark Matter and Energy 

5.1 Degrees of Freedom 

A dark component interacts with ordinary matter only through gravity and hence its observable 
properties are completely specifled by the degrees of freedom in its stress-energy tensor. We have 
seen that without loss of generality these can be taken as the elements of the symmetric 3x3 stress 
tensor. Two of these elements represent the scalar degrees of freedom that influence the formation 
of structure. In the background only one of these remain due to isotropy leaving only the pressure 
Pd oi the dark component or equivalently the equation of state wd = Pd/pd- Let us generalize 
the concept of equations of state to the fluctuations in the stress tensor. 

5.2 Generalized Equations of State 

It is convenient to separate out the non-adiabatic stress or entropy contribution 

PD^D = SpD - c\,JpD , (103) 

where the adiabatic sound speed is 



„2 ^ 



PD 1 WD (a\ 

WD - TTT— - . (104) 



PD ?,l + WD\a 



Therefore, pD = wdPd does not imply 6pD = wd^pd and furthermore if Td ^ 0, the function 
WD{a) does not completely specify the pressure fluctuation. For the dark energy, wd < and is 
slowly- varying compared with the expansion rate (a/a) such that c|p^ < The adiabatic pressure 
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fluctuation produces accelerated collapse rather than support for the density perturbation. There- 
fore a dynamical dark energy component must have substantial non-adiabatic stress fluctuations 
7^ to be phenomenologically viable. 

One cannot simply parameterize the pressure fluctuation by a non-adiabatic sound speed 
Spo/SpD since this is not a gauge invariant quantity and the non-adiabatic stress fluctuation is. 
The gauge covariant approach allows us to define the equation of state covariantly. 

As we have seen the comoving gauge gives the closest general relativistic analogue to non- 
relativistic physics. The generalization of the comoving gauge to an individual dark component has 
the conditions 

B = VD, Ht = 0, (105) 



which we will call the rest gauge since Tf = 0. Taking 

cl , (106) 



Sp 

Tp 



rest 



we obtain via gauge transformation 



OPD = Cr,dpD - {cdPd - Pd) — — 

= clSpD + 3(1 + WD)%Dicl - claf-^^ , (107) 

yielding a manifestly gauge invariant non-adiabatic stress 

(108) 



d V£) — B 
OPD + 3(1 + WD)-PD 



o k 



The anisotropic stress can also affect the density perturbations. A familiar example is that 
of a fluid, where it represents viscosity and damps density perturbations. More generally, the 
anisotropic stress component is the amplitude of a 3-tensor that is linear in the perturbation. A 
natural choice for its source is kvr,-, the amplitude of the velocity shear tensor However it must 

also be gauge invariant and generated by the corresponding shear term in the metric fluctuation 
Rt- Gauge transforming from a gauge with E.t = yields an invariant source of kvo — Ht- These 
requirements are satisfled by 

^ + Z^WDT^D = '^c]jAkVD-HT)- (109) 



5.3 Examples 

Cold dark matter provides a trivial example of a dark component. Here stress fluctuations are 
negligible compared with the rest energy density and so wcdm = = c|)j^ = 0. Scalar field 
dark energy provides another simple example. Here the competition between kinetic and potential 
energy can drive < 0. As we have already seen a slowly rolling scalar field has = 1 by virtue 
of the absence of scalar field fiuctuations in the rest gauge. In this case the fluctuations bear only 
the kinetic contributions and so the relationship is exact. 

The non-adiabatic sound speed c|) is also useful in characterizing A;-essence, a scalar field with 
a non-canonical kinetic term. Since in the rest gauge the potential fluctuation vanishes, the sound 
speed directly reflects the modification to the kinetic term. A special case occurs when the kinetic 
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Type 




r2 
u 


r2 

^ Tin, 

uv 


CDM 











A 


-1 






Massless v 


1/3 

/ 


1/3 


1/3 


Tight coupled 7 


1/3 


1/3 





Hot /Warm DM 


1/3^0 


same 


same 


Quintessence 


variable 


1 





k-essence 


variable 


variable 





Phantom energy 


< -1 


1 





Decaying v 


1/3^0^ 1/3 


same 


same 


Axions 





small 





Fuzzy Dark Matter 





scale dependent 






term has the wrong sign. Here < — 1 and the energy density increases as the universe expands. 
This type of matter has been dubbed phantom dark energy. 

The viscosity term is relevant for dark radiation components. Massless neutrinos can be ap- 
proximated by Wi, = cl = = 1/3. A massive neutrino has an equation of state that goes from 
Wmu = 1/3 to as the neutrinos become non-relativistic. Fitting to the numerical integration of 
the distribution gives 

Wm. = l[l + {a/a^rfP]-'^\ (110) 

with p = 0.872 and Onr = 6.32 x 10^^ /Q^h'^ . We can model its behavior as c^^ = c^,^^ = wd- 

A summary of the phenomenological parameterization of various particle candidates are given 
in the table. 



5.4 Initial Conditions 

With the equations of motion for the dark components defined, all that remains is to specify the 
initial conditions. Conservation of the comoving curvature outside the Hubble length allows us to 
ignore the microphysics of the intermediate reheating phase between inflation and the radiation 
dominated universe. One can then simply take the initial conditions for structure formation to be 
in the radiation dominated universe. Accounting for the neutrino anisotropic stress as described 
in the previous section, the initial conditions for the total perturbations in the comoving gauge 
become 



5 = Asikv)\, 
U = Au{kri)H, 



(111) 
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where 5 = Sp/p, and r] is now understood as the conformal time elapsed after the inflationary 
epoch. The constants are 

41 + 2/./15 

91 + 4/,,/15^ ' 



1 1 



3 1 + 4/,/15 ' 



4 1 

M = (112) 

15 1 + 4/^/15 ' ^ ' 



and 



fu = ^^. (113) 

P7 + Pu 

The dark component initial conditions can then be determined by detailed balance 

6j = Asjikr^fC 
vj - V = AAvj{kvfC, 

Uj = AnAkvfC (114) 



with 



r2 



Auj = ^Au , 

WJ 

, ^[2 - 3(^, - 4)](/. - - [2(4 - i) + (^., - 4)](1 + A/,) ^ 

(4 - 34) [2 - 3(^,-4)] + 9(4- 4 Jc^, 'a:^^"'^ 



These initial conditions apply to all particle species J, including the photons and neutrinos, save 
that for the baryons rapid Thomson scattering with the photons sets Vb = v-y. In fact for most 
numerical purposes one can simply set vj — v = in the initial conditions and let the velocity 
differences arise dynamically in the integration. The initial conditions in an arbitrary gauge can be 
established from these relations and the gauge transformation properties of the perturbations. 



6 Cosmic Microwave Background 
6.1 Boltzmann Equation 

The gauge covariant formalism is useful for CMB studies as well. The interpretation of CMB 
anisotropy formation is simplest in the Newtonian gauge where the manifestations of gravitational 
redshift and infall correspond to Newtonian intuition. The numerical solution of the perturbation 
equations is best handled under the comoving or synchronous gauge where the fundamental pertur- 
bation variables are stable. The link to inflationary initial conditions is best seen in the comoving 
gauge. The gauge covariant approach allows one to calculate in one gauge and interpret in or relate 
to another. 

Most of the physics of CMB perturbation evolution is contained in the general discussion of ^ 
However the cosmic microwave background differs from the dark components in that it undergoes 



27 



interactions with the baryons that can exchange energy and momentum. Thus the conservation 
law for its stress energy tensor must be supplemented with an interaction term. CMB observable 
properties are also not the energy density and velocity perturbations but the higher order angular 
distribution of its temperature and polarization. 

Generally the CMB is described by the phase space distribution function photons in each of the 
two polarization states /a(x, q, rj), where x is the comoving position and q is the photon momentum. 
The evolution of the distribution function under gravity and collisions is governed by the Boltzmann 
equation 

-^fa,b = C[fa,h], (116) 

where C denotes the collision term. 

In absence of collisions, the Boltzmann equation becomes the Liouville equation. Rewriting the 
variables in terms of the photon propagation direction n, 

^/„(x, n, q, ri)=fa + nWifa + q^fa = • (117) 

The last term represents the gravitational redshifting (or Sachs- Wolfe effect) of the photons under 
the metric and is given by the geodesic equation as 

^ = - - - In'n^HTH -Hl + n'Bi - n'ViA . (118) 
q a 2 ■' 

We have in fact already derived these gravitational effects in the general covariant perturbation 
formalism. To establish this fact note that the stress energy tensor of the photons is the integral 
of the photon distribution function over momentum states 

j^ni^ I d^q q^q 



The Liouville equation then expresses the conservation of the stress-energy tensor. Given that the 
CMB distribution is observed to be close to blackbody, it is suffices to calculate the evolution of 
these integrated quantities. In particular let us define the temperature perturbation 

e(x, n, r?) = \6, = ±-J^{f, + f,)-l (120) 

and likewise for the linear polarization states Q and U as the temperature differences between 
the polarization states. The redshift terms then become the metric terms in the continuity and 
Navier-Stokes equations. 

6.2 Eigenmodes 

As in the case of the purely spatial perturbations, we decompose the temperature and polarization 
in normal modes of the spatial and angular distributions. For the A;-th spatial eigenmode 



e(x,n,ry) = ^e('")Gr(k,x,n), 
[Q±iU]{^,n,v)=Y^[Et^±iBi"'^]±2Gf{K^,n). (121) 
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They are generated by recursion {qG = G) 



21 + 



qms 



where 



(122) 



q"^ = k^ + {\m\+l)K . 



(£2_^2)(£2_g2) 



£2 ■ 



(123) 



The lowest order modes begin the recursion and are related to the spatial harmonics as 

qCtj- -n ...n I iC^ij...i|^| , 
±2^^ oc (mi ± im2f{7ni ± im^TQ^^, , 



(124) 



where liii , 1112 span the plane perpendicular to fi. The normalization is set so that 

,GT{0,0,n) = i-i) 



(125) 



Here the spin spherical harmonics ^Yf^ are the eigenfunctions of the 2D Laplace operator on a 
rank s tensor. They are given explicitly by rotation matrices as 



2i + l 
47r 



(126) 



The meaning of these modes becomes clear in a spatially flat cosmology. Here the modes are 
simply the direct product of plane waves and spin-spherical harmonics 



Air 



Gr(k,x,n) = (-^)Y y™(fl)exp(^k•x 



±2G,-(k,x,n) = (-i) 



2e + 

2£+l 



±2yr(n)exp(ik-x). 



(127) 



The main content of the Liouville equation is purely geometrical and describes the projection of 
inhomogeneities into anisotropies. Photon propagation takes gradients in the spatial distribution 
and converts them to anisotropy as 



47r, 



fi . Ve^""-^ = in ■ ke*^-^ = -jkY^{n)e^^-'' . 
This dipole term adds to angular dependence through the addition of angular momentum 



(128) 




yOyrn 



v/(2^ + l)(2^-l) 



K 



v/(2^ + l)(2^ + 3) 



V 



e+i ' 



(129) 



where kJ^ = = y/P — 171"^ is related to the Clebsch-Gordon coefficients. 
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6.3 Collision Term 

The dominant collision process for CMB photons is Thomson scattering off of free electrons which 
has the differential cross section 

^ = A|E'.E|V, (130) 

where E' and E denote the incoming and outgoing directions of the electric field or polarization 
vector. 

To evaluate the collision term we begin in the electron rest frame and in a coordinate system 
fixed by the scattering plane, spanned by incoming and outgoing directional vectors — n' • n = cos P, 
where (3 is the scattering angle. Denoting as the in-plane polarization temperature fluctuation 
and @± as the perpendicular polarization state, we obtain the geometrical content of the transfer 
equation 

Gy oc cos^ f3 ej| , oc e'^ , (131) 

where the proportionality reflects the scattering rate 

r = n^aTa . (132) 

To calculate the Stokes parameters in this basis, we also need to calculate the polarization states 
with axes rotated by 45° 

El = (E|| + E^) , E2 = (E|| - E^) , (133) 

yielding the transfer 

61 (X |Ei • Eil^e'i + |Ei • Ease's 

oc ^(cos/3 + 1)2g; + h^cos 13- if Q'2 
82 (X |E2- Ease's + |E2-Ei|2e; 

cx ^(cos/? + + ^(cos (3 - ife'^ . (134) 

Now the transfer properties of the Stokes parameters 

6^^(011+0^), Q ^ 1(011 -e^), C/^ 1(61 -62) (135) 

arranged in a vector T = (0, Q + ill, Q — ill) becomes 

TcxS(/3)T', (136) 

„ , cos^/? + l -isin^/3 -isin^/3 
S(/3) = -| -isin^Z? i(cos/3 + l)2 i(cos/3-l)2 



-^sm^(3 i(cos/3-l)2 i(cos/3 + l)2 



4 

where the normalization factor is set by photon conservation in the scattering. 



30 



Perturbation Theory 



Finally convert the the polarization quantities referenced to the scattering basis to a fixed basis 
on the sky by noting that under a rotation T' = R('(/')T where 

1 

R(V;) = I I , (137) 



e'^'^ 



giving the scattering matrix 
R(-7)S(/3)R(a) = 



r2°(/3, a) + 2V5roO(/3, a) -^lY~^fi, a) -^jlY^iP, a) 

-x/6_2K0(/?,«)e-'*^ 3_2y2-^(/3,a)e-2,7 3 _2K2(/3, a)e-2»^ 



where a, 7 are the angles required to rotate into and out of the scattering frame. 
Finally, by employing the addition theorem for spin spherical harmonics 



^ ,,yr*(n') sjrii^) = (-l)'^"'^y^^ ,^r,-^^(/3,a)e-^^ (139) 

m 

the scattering in the electron rest frame into the Stokes states becomes 
a„[T]=f j ^R(~7)S(/3)R(a)T(n') 

2 

= f / ^(e',0,0) + ^f / dn' ^ p('")(n,n')T(n'), 

m— — 2 

where the quadrupole couphng term is 



(140) 



Y^»{n')Y^{n) ~^l^Y^*{h')Y^{n) -.^'^ _^Yi^*{h')Y^^{n) 
P(™)(h,n')=| -V61T*(h')2>^2"(A) i^Yr{n')^Y^{n) 3 _^Yr*{n') ,Yrin) ]■ (141) 
-V6Yr*{n')_^Yr{n) 3 ^Y^* (n') _^Yr (n) 3 _2Yr*{n') _^Yr{n) 

The full scattering matrix involves difference of scattering into and out of state 

C[T] =an[T]-Cout[T]. (142) 

In the electron rest frame 

C[T] = f f — (e', 0, 0) - tT + Cp[T] (143) 
J 47r 

which describes isotropization in the rest frame. All moments have e^'^ suppression except for isotropic 
temperature Oq. Here Cp is the P or £ = 2 term of Eqn. 1140(1 . 

Transformation into the background frame simply induces a dipole term 

C[T] = f ("n • Vfa + f ^e',0,o \ -tT + Cp[T], (144) 



yielding the final form of the collision term. 
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6.4 Temperature-Polarization Hiercirchy 

The Boltzmann equation in normal modes then becomes 



E^' = k 



B^r'^ = k 



21- 

2 4 



(m) 
«-l " 

(m) 



2^ + 3^+1 



2m „(m) _ 2K^i 



2^ + 3 



2^-1 £(^+1) 
where the gravitational and scattering sources are 



2m ^(m) _ infj^x 



2i + 3 



(145) 





V 

-7-\/6P('")5^,2 , 



f _ Jj(±2) 



with 



(146) 



(147) 



The physical content of the coupling hierarchy is that an inhomogeneity in the temperature or polariza- 
tion distribution will eventually become a high multipole order anisotropy by "free streaming" or simple 
projection. 

6.5 Integral Solution 

Since the hierarchy equations simply represents geometric projection, their implicit solution can be written 
as the projection of the gravitational and scattering sources at a distance. This operation proceeds by writing 
the normal modes themselves in spherical coordinates, 

e 

±^GZ = ^(-i)V47r(2^+l)[45 ± J](fc, D) ^,Yr{n) , 



(148) 



where D = rjo — r]. Summing over the sources 



2^+1 

El'"\k,m) 

21+1 
Bt\k,vo) 

2e + i 



/ d^e-^4")4:;(fc,D), 

° is 

/ ^ / N 

/ d^e-^£:i:)/3ff(fc,i3). 

^ 



(149) 



Note that the polarization has only an = 2 source. 

In a flat cosmology, the radial projection kernels are related to spherical Bessel functions 



gik.x ^ ^(_,)^ V4^(27TT)j^(fcZ?)r/(n) 



(150) 
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by the recoupling of the "spin" angular dependence of the source £s to the "orbital" angular dependence of 
the plane waves. For example 

a'^^^{k,D)^jfikD), 
af^{k,D)=f,{kD), 



^'^'^^-y 8(£-2)! (kD)^ ' 
f3if{k,D)=0. (151) 

In a curved geometry, the radial projection kernels are related to the ultraspherical Bessel functions by the 
same coupling of angular momenta. 

6.6 Power Spectra 

The two point statistics of the temperature and polarization fields are described by their power spectra 



,_2 fdk^lf{X^;^^2^ n.o^ 



where X,X' G Q,E,B. 

7 Epilogue 

How to conclude Lecture notes in which no conclusions have been drawn? I leave you instead with a thought: 

What goes on being hateful about analysis is that it implies that the analyzed is a completed 
set. The reason why completion goes on being hateful is that it implies everything can be a 
completed set. 

-Chuang-tzu, 23 
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